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SUMMARY 


The sound -pres sure field of a rotating propeller in forward flight 
in free space is analyzed by replacing the normal-pressure distribution 
over the propeller associated with thrust and torque by a distribution 
of acoustic pressure doublets acting at the propeller disk and subject 
to uniform rectilinear motion. The basic element used to synthesize the 
field is the pressure field of a concentrated force moving uniformly at 
subsonic speeds, for which an expression generalizing one of Horace Lamb's 
for the fixed concentrated force is given. This result is presented both 
for the moving and for the fixed observer. The strength of the doublet 
distribution is related to the thrust and torque distribution in a con- 
venient but approximate way. The sound field is expressed by integration 
over the propeller disk, and also by integration over an effective ring, 
and is given both for the near pressure field and, in a simpler form, 
for the far field. Known results for the zero-forward-speed case present 
themselves in the special case of Mach number M = 0. Some illustrative 
examples are calculated and discussed. 


INTRODUCTION 


The rotating propeller is the source of an intense sound-pressure 
field which can be associated with the periodic reactions on the medium 
arising from the distribution of pressure rotating along with, the blades. 
This pressure distribution consists in part of a distribution due to 
thickness of the blades, whose resultant force in subsonic potential flow 
is zero, and in part a distribution due to angle of attack and camber 
of the blades whose integrated effect includes the induced drag and corre- 
sponds almost wholly to the thrust and torque distribution over the blade. 
Another source of propeller noise may be associated with flow separation 
and with friction or shear due to the boundary layer j both effects lead 
to vorticity shed into the wake and hence the designation vortex noise. 

The vortex noise and the noise due to thickness (where wave drag is not 
a large factor) are, however, for actual propellers normally of a consid- 
erably smaller magnitude than the rotational sound due to torque and 
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thrust; hence only the latter effect will be considered in the present 
work. 


A large number of investigators have studied various phases of the 
determination of the sound or noise field of rotating propellers. In 
addition to the references cited, a bibliography is included of repre- 
sentative work on this subject. A simplification that has frequently 
been made in propeller -noise investigations is to limit the consider- 
ations to the static or standing propeller. The work of Gutin in 1936 
(ref. l) represents a development of this type that makes possible a 
satisfactory prediction of the amplitude of sound pressure due to thrust 
and torque of a propeller rotating on a stand in still air. Although 
Gutin' s theory is applicable to the near oscillating-pressure field of 
the propeller, his results in reference 1 are limited to the determi- 
nation of the fundamental and the first few harmonics at a distance far 
from the propeller, that is, several diameters away. The determination 
of the near pressure field, however, has been of concern both from struc- 
tural and physiological considerations. Hubbard and Regier (ref. 2) 
extended the application of Gutin' s work to .describe the oscillating 
pressure field and to determine the amplitude of noise at points near 
the propeller, in some cases within a blade chord length from the tip. 
They investigated analytically the effect of several of the parameters 
that enter in the theory and also gave comparisons with experiment which 
were quite satisfactory. 

The existing theoretical work has found useful application for 
static conditions and for conditions of low forward speed, for example, 
near take-off. For conditions of high forward speeds, however, many 
pertinent questions have arisen as to the possible effects of the forward 
speed on the oscillating pressure field of the propeller. A few investi- 
gators have examined phases of this problem; Bryan, Hart, Shirokov, 
Blokhintzev, and particularly Kuessner and Billings (see bibliography) 
may be mentioned, but most of this work seems incomplete or difficult to 
apply. It appeared desirable, especially for applications, to reexamine 
the theoretical problem for the propeller at forward speed so as to have 
it arise as a straightforward generalization of existing work for the 
zero -forward-speed case. The purpose of the present paper is therefore 
to extend the theory of Gutin (ref. l) and the work of Hubbard and 
Regier (ref. 2) so as to include the effect of subsonic forward speed 
of the propeller on the near and the far oscillating sound -pressure 
field caused by torque and thrust of the propeller. 

This paper includes the following material: (a) The sound- pressure 

field associated with a uniformly moving concentrated force is given with 
details of the explicit development contained in the appendix. This 
result serves as a basic element for synthesizing the total pressure 
field. It is expressed both for the case of the observer or field point 
considered to be moving along rectilinear ly and uniformly with the pro- 
peller in free space and for the case of the observer considered fixed 
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and the propeller in uniform flight, (b) The disturbance forces associ- 
ated with the pressure distribution acting on the medium in the plane of 
the propeller disk are presented in the manner similar to that of Gutin. 
(c) The sound -pressure field resulting from the combined use of (a) and 
(b) is shown in the form requiring integrations over the propeller disk 
and in the simpler form making use of an effective propeller radius. 
Approximations valid for the far field, which reduce to Gutin 's results 
for the case of zero forward speed, are also given. (d) Some numerical 
examples are calculated and discussed. (e) Remarks on the usef uln ess 
and limitations of the analysis. 


SYMBOIS 

A(r) chordwise distribution of thrust acting on a radial element 

of a propeller blade 

B number of propeller blades 

b(r) width of propeller blades 

c velocity of sound 

Cp power coefficient, p/pn^D^ 

Crp thrust coefficient, T /pn^D^ 

D propeller diameter 

F ( r ) chordwise distribution of forces perpendicular to the thrust 

of a propeller blade and giving rise to the torque 

F x> F y> F z components of force vector F 

Bessel function of first kind with index raB 

Mach number, V/c 
order of harmonic 
propeller rotational speed, rps 


k-2 

c 


J mB 


M 

m 


n 
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P 
P 

P 

P Q 
P T 

^rms 


power 

pressure 

pressure magnitude 
pressure due to torque 

pressure due to thrust 
root-mean-square pressure 


Q torque 

R length of propeller blades 

R e effective length of propeller blades 

r radius to a blade element 

r, 0 polar coordinates in yz-plane 

s = f - x-l ) 2 + P 2 [(y - yp) 2 + (z - .if] 

S Q = ^/x 2 + p2y 2 


s = ]f(x - x x ) 2 + (y - y x ) 2 + (z - z 
s Q = \/x 2 + y 2 

Sq, polar coordinates in xy-plane 


T 

t 

V 


thrust 

time 

forward velocity 


x, y, z 


Cartesian coordinates 
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3 = 


£ 


M 2 


Pr 


blade angle 


q propeller efficiency 

p fluid density 


M(x - x^) + S 



T 


_b_ 

rn 


T 


O 


ft 


CD 

CDi 


period, 2rt/Bft 

angular velocity 

frequency of mth harmonic, mBQ 

fundamental frequency. Bn 


ANALYSIS 

Extension of a Formula in Lamb's "Hydrodynamics" for the 
Sound-Pressure Field of a Fixed Concentrated Force 
to That of a Moving Force 


On the basis of acoustic considerations of the classical hydrodynamic 
equations. Lamb (ref. 3 ) gives the pressure at any field point x, y, z 

associated with an external periodic force Ae icut acting in the x-direc- - 
tion and concentrated at the location x^, y-^, z-^ as 


a a 

8x s 


( 1 ) 


where 

= )j(x - x-l ) 2 + (y - yj 2 + (z - Zl ) 2 


s 


( 2 ) 
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The concentrated force may be regarded as stemming from the limit of a 
distribution of an increasing pressure difference over a decreasing 
area, whose product in the limit of zero area is equal to the force. 

For an arbitrary concentrated force of components F x , F y , F z the result 

generalizes as Iamb indicates to 


_1_ 

a F x 

(-■§) 

a F y 

4 _ - ■ 

(-1) 

a F z 

(*-i) 


4rt 

dx 

s 


s 

dz 

s 


(3) 


This formula shows that the pressure field obtained from the concentrated 

force has the character of a doublet or a dipole 1 ) in its periodic form 
it has played a central role in the development of the noise and pressure 
fields of propellers due to the- propeller torque and thrust. However, 
the formula refers to a concentrated force and to coordinates fixed in 
space. To treat the case ihvolving effects of forward speed, it is expe- 
dient to obtain the required extension of the preceding formula. 

The extension of equation ( 3 ) required for the case of a uniformly 
moving concentrated force is given in equation (4), which follows. 

Details for its derivation are supplied in the appendix. let the con- 
centrated force have components F x , F y , F z ) let it be moving uniformly 

with velocity V in the direction of the positive x-axis; and let the 
coordinate system also move uniformly with the same velocity. Then, 
with x, y, z now denoting the field point referred to this coordinate 
system and with x^, y^, z-j_ denoting the coordinates of the concentrated 

moving force 



a F y(* - 1) , a *•*(* - 1) 

by S bz S 


(4) 


where 


S = )j(x - x-l ) 2 + 3 2 [(y - y±) 2 + (z ~ Zi f\ 
M(x - x-l) + S 


3 = VI - M 2 


(5) 


1 Since the acoustical pressure is p = -p the dipole is an accel- 
eration doublet or the time rate of change of a flow doublet. 
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Figure 1 illustrates the geometric meaning of the quantities s, S, 
and cr. The force located at C^x-^y-^z-jJ and the field point at 

Q(x,y,z) are both in uniform motion with velocity V in the x-direc- 
tion. The distance OQ = s. The influence at Q at time t stems 
from the action of the force when it was at position P, where P is 
obtained from the relation OP/OQ = M, the forward-speed Mach number. 

The distance PQ gives then directly the phase_radius a. The perpen- 
dicular dropped from 0 onto PQ determines QR, which is equal to the 
amplitude radius S. For M = 0 both S and o reduce to the ordinary 
radius s. Comparison of equations (3) and (4) shows that the effect of 
the forward speed leads to replacement of s by S in the amplitude and 
by a in the phase. 

The force variation of main concern herein is that of a harmonically 
periodic force F having components V F z varying in time as e iaj t. 

Equation (4 ) becomes 


P 





a\ e - ikg 

dz/ S 


(6) 


where k = cjd/c. Equation (6) is the basic relation to be ma de use of in 
the subsequent analysis. 


Disturbance Forces in the Propeller Plane 

Consider as does Gutin (ref. l) a propeller disk oriented so that 
the axis passing through the center of the disk coincides with the x-axis 
and let the propeller be assumed to rotate, as in figure 2, in the 
yz -plane (x = 0) with positive values of x corresponding to points 
ahead of (and negative values of x corresponding to points behind) the 
propeller disk. The propeller is considered to move uniformly with 
velocity V in the positive x-direction. Points in the propeller disk 
are designated by (0,y^,z.j^, or* in polar coordinates (as shown in 

fig. 2), by 


y 1 = r cos 9 
z^ = r sin 9 

For definiteness, let the propeller be rotating counterclockwise as seen 
by an observer looking into the propeller toward the slipstream. 
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Each element of the propeller Is acted on by the surface pressure 
distribution and this distribution may be resolved into a thrust force 
in the direction of the axis of rotation (the x-axis) and into a force 
associated with the torque which acts about the axis of rotation opposing 
the rotation. Equal and opposite reaction forces to these are exerted on 
the medium. The points of application of these forces are imagined to 
act in a single plane designated as the plane of rotation or as the pro- 
peller disk. 

The reaction of the surface pressure distribution of the rotating 
propeller on the medium at any instant is to be replaced by fixed 
periodic forces acting at the- propeller disk. The propeller disk itself 
will be considered covered with the necessary singularities in the pres- 
sure or acoustic radiators of proper strength and harmonic content to 
correspond to this normal-pressure distribution. These singularities 
will be seen to be acceleration sources (acoustic radiator of zero order) 
for those forces which act symmetrically on both sides of the blade, that 
is, whose net force over each element of the blade is zero, and to be 
acceleration doublets (acoustic radiator of order one) for those forces 
acting antisymmetrically, that is, whose net force over each element of 
the blade corresponds to the difference in pressure over both sides of 
the blade. The doublet distribution is that which is needed to repre- 
sent the thrust and torque distribution, in particular that part of the 
thrust and torque distribution associated with pressures acting normal 
to the blade surface. This part is practically all of it, except that 
arising from skin friction. (The effect of blade thickness may be taken 
into account by introducing either flow sources or acceleration sources: 
the flow sources would lead to the sound-pressure field due to thickness 
by the classical "piston" effect of the moving blade on the flow field; 
the acceleration sources would deal directly with the assumed or known 
contribution to the pressure distribution over the blade due to thickness.) 

Consider an element of the propeller at distance r from the axis; 
let dr be its radial length and b its width measured in the projec- 
tion onto the plane of rotation. Let the forces acting on the propeller 
element on each blade be A(r )dr in the axial flight direction and 
F(r)dr in the direction opposed to the direction of rotation. Equal 
reaction forces acting opposite to the flight direction and in the direc- 
tion of rotation, respectively, are exerted on the medium. These quanti- 
ties are related to thrust T and torque Q by the relations 


^ = BA(r) 
dr 

= BrF(r ) 


(8a) 
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or 



(8b) 


where B is the number of blades and R the length of each blade. 

The periodic impulse or reaction experienced at any element of the 
disk may be expanded in a Fourier series. For simplicity let the element 
considered first be located on the radial line 9 = 0 (the y-axis ) and 
afterwards be located arbitrarily. To start with, assume that the forces 
are uniformly distributed over the projected width b of the blade 
element, that is, that the distribution of pressure difference over the 
blade chord is rectangular. Then in the area element r dr d9 of the 

plane of rotation, the forces A(r)dr r ^ and F(r)dLr r ■ act on the 

b b 

medium during the time interval in which this element is eclipsed by the 
projection of the propeller element. If the overlapping of the element 

starts at t = 0, it will end at t = r = — and the overlapping of the 

rft 

element will start again by the next blade at t = t q = — , where Q is 

bh 

the angular velocity of the propeller. The rectangular -type forces 
experienced at the element of the disk located at 9 = 0 by its periodic 
eclipse by the blade may be developed in a Fourier series: 


and similarly 


Fx*(t) = 


F x *(t) = 


F 2 *(t) = 


A(r) — dr for 0 < t < r 

JO for t < t < r 0 

00 

Aq + Aj n cos(mBftt - e m ) 

m=l 



r dQ 
b 


dr 


for 0 < t < t 


[0 for t < t < t q 

00 

F 2 *(t) = B 0 + ®m cos(mBnt - T) m ) 
m=l 


( 9 ) 


( 10 ) 
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where 


Bm 


2 r y 

— A(r) — sin rrur — dr d0 

nut b t 0 

— F(r) — sin rmt — dr d0 

rrur "b T n ! 


t mB^T 

% = "w — = — 

T o ^ T o 


( 11 ) 


The constants A q and Bq which correspond to the instantaneous average 

thrust and torque over the blade element and to the associated momentum 
shed into the slipstream do not give rise to sound and need not be 
expressed. The phase angles e m and q m , which are small for the lower 

harmonics, are needed to preserve generality of discussion since Ajjj cos e m 
corresponds to the coefficients of pure cosine terms and sin e m to 

coefficients of pure sine terms in the Fourier series. For the assumed 
rectangular distribution, e m and may be observed to be equal to 

zero if the origin of time is chosen at the overlapping of the center 
line of the blade by the y-axis since then only a pure cosine series 
suffices to express the distribution. 


The general formulation of interest, where the disk element r dr d9 
is located at an angle 0, may be expressed directly from these results, 
for, on a second area element r dr d9 shifted with respect to the first 
by the angle 0 in the rotational direction, there act periodic forces of 
the same magnitude, but retarded by the time 9 /n . The corresponding 
Fourier developments are 


F x (t) = F ± *(t 
F 2 (t) = F 2 *(t 



( 12 ) 


The quantity 


mjtT _ mBfir _ mBb 
Tq 2 " 2r 


( 15 ) 
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is small for the lower harmonics, especially for blade elements not near 
the hub. Blade elements near the hub, where r is small, are eliminated 
from consideration since they contribute very little to the air forces. 

One may then replace in equation (ll) sin by 22LL so that, noting 

To T n 

the relations in equations (8a), 


An ” f A(r)dr 

%“f F (r)dr 


This approximation, which weights 

monies equally, may be observed to correspond to the case of the thrust 
and torque distributed over a zero blade width (that is, to the mathe- 
matical pulse sometimes termed the Dirac delta function). This approxi- 
mation becomes relatively less valid when exceeds about 

t 0 4 

sin = 0.707» 7 — = O. 785 ) or when 4m — exceeds about unity; that is, 

4 4/ t q 

2 . 0 , 

when the order of the harmonic exceeds — or, roughly, l/solidity. 

Bb 

Moreover, the assumption of a uniform rectangular distribution of 
the forces across the blade has been made for convenience, and other 
distributions may be treated if desired. Some of the possible errors in 
the assumption leading to equation (l4) are discussed briefly by Gutin 
in reference 1 and are shown to be generally negligible for the lower 
harmonics. Regier and Hubbard (ref. 4) also discuss this assumption in 
an illuminating manner. Figure 3> which is taken from reference 4, 
shows the relative harmonic content of different assumed distributions 
of the same total load for: (a) the mathematical sharp pulse, (b) a 

"b 

triangular "hat" pulse = 0.03, (c) a rectangular pulse = 0.03, 

2itr 2nr 

and, (d) a rectangular pulse — = 0. 06. Blade widths of actual pro- 

2nr 

pellers tend to be between cases (c) and (d) for the most effective parts 
of the propeller disk. The sha.rp-pulse assumption generally tends to 
overestimate the magnitude of the higher harmonics. When the proper 
distribution is known, appropriate correction factors for the required 
harmonic may, of course, be applied to the magnitude of the results given 
by the pulse solution. 


d9 


d9 


dT 

— dr d0 
dr 


- — dr d0 
r dr 


(14) 


the harmonic content of all the har-- 
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The formulation in equation (12) need not be limited to. the one in 
which the propeller force distribution is uniform throughout the pro- 
peller cycle. If interest should be attached to a nonuniform distribu- 
tion, as occurs for the propeller yawed or pitched with respect to the 
■ flight path or for one experiencing interference effects, it is readily 
possible to allow for these effects by permitting the amplitude of the 
distribution to become a function of 0. 


The Sound-Pressure Field 

The general expressions for the resolution of the forces associated 
with the thrust and torque on a radial blade element (for example, 
eqs. (12)) may be put in the usual convenient complex form (whose real 
part may correspond to the formulation of interest): 


F 1 - Ao + zL 

m=l 

00 

F 2 = B o + X! 

m=l 


V 


i(a)t-mB9-e ni ) 


> 


(15) 


®m e 


i(u)t-mB9-i} m ) 


where go = mBW and where for the special case of the rectangular pulse 
distribution the coefficients are given as in equation (11 ). 

With the use of the approximations for and given in equa- 

tion (14), the periodic forces acting on the medium at the element r dr d9 
of the disk can be expressed for any given harmonic in the forms 


F 

r x,m 

F 

F 

x z,m 


1 dT i(oot-mB9 ), JA 

= e v dr d9 

n dr 

1 1 dQ . . i(o>t-mB0 ) JA 

= "t— — sin 0 e v dr d0 

it r dr 

1 1 dQ . i(cut-mBe), JA 

= — ■ - cos 0. e v ,dr d9 

it r dr 


> 


(16) 


In the derivation of equation (l6), the phase angles e m and T] m , which 

can generally be neglected for the lower harmonics, have been put equal 
to zero, a value which corresponds in the case considered (as has been 
mentioned) to the choice of 0=0 to correspond to the overlapping of 
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the center line of one of the blades at the time ' t = 0. The index m 
indicates the order of the harmonic considered. As discussed in the 
preceding section, appropriate factors depending on the harmonic num- 
ber m, functions of the radius, and, for unsymmetrical loading, func- 
tions also of 9 may (if known) be applied to the terms in equation (l 6 ). 

Pressure relations involving integrations over the, propeller disk . - 
The pressure at any point in the free-space field produced by these com- 
ponents of the uniformly moving periodic forces F x , F y , F z is given in 

equation ( 6 ). With the use of the components of force given in equa- 
tion (l 6 ) and by integration over the propeller disk, the total oscil- 
lating pressure p for any given harmonic m (the index m will be 
hereinafter dropped) is obtained from equation ( 6 ) as a sum of the pres- 
sure due to thrust p T and that due to torque Pq as follows: 

P = P T + P Q (17) 


where 


P 


T 


1 f R r 2n dT e i(ojt-mB0) 

bit 2 J 0 J 0 dr 


A e ~ ikg 

cbc S 


dr d 0 


(18) 



1 r R r 2n 1 dQ 

4rt 2 J 0 J 0 r & 


e i(u>t-mB9) 1 8 e~^ kcT 
r 89 S 


dr d9 


(19) 


and where, in the expression for p n , the operation A A has been used 

W r o9 

for convenience to replace the equivalent operation sin 9 A - cos 9 A. 

oy 8 z 

The indicated differentiations with respect to x and 9 can be carried 
out by use of the expressions for a and S given in equations ( 5 ) 


8 e"^ kg _ e ^ kq / ikM ikx x \ 

^ S S V 3 2 " ’ s2/ 


1 _ 8 _ e ~ ikCT = e - lk(T 

r d9 S ~ S 


ik 

S 





sin 9 - z cos 9) 


( 20 ) 
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Hence, 


3 = r R r 2n dT e -imB0 £^/ikM + ikx + xJW d0 (21 ) 

' T kit 2 J 0 J 0 ^ s \p 2 (B 2 S s 2 / 


and 


e tot f E r 2 ” 1 aq -UnB9 el k0 

PQ ' 4* 2 J 0 J 0 r dr 3 


(f*S) 


(y sin 0 - z cos 0 )dr d0 


(22) 


Another somewhat simpler expression for p^ can be obtained from equa- 
tion (19) by integration by parts with respect to 0 


P « 4« 2 J 0 Jo rdr s 


(23) 


Pressure relations involving effective ring approximation . - Appreci- 
able simplification may be achieved for calculation purposes by making 
use of the approximations inherent in the assumption of an effective pro- 
peller radius R e so that the integration with respect to r is avoided. 

Equations (21) and (23) then reduce to 


.icot p2rt 


Pqi - 


4it £ 


r T /ikM + ikx + xV 

Jo \p 2 P 2 S s 2 / 


i -i(mB0+ka) 


d0 


and 


acut p2n 


P Q , 2 


Q -i(mB0+ka ) 

J±_ imp d0 


4rt* R e 


(24) 


(25) 


where, in both S and 0, the points y^ and have the values 

y^ = Re cos Z 1 = R e sin 0 * In e ^ ec " t: > equations (24) and (25) imply 

that the propeller disk has been replaced by an annular ring in which 
the entire thrust and torque are concentrated. The effective radius of 
this ring varies somewhat, as calculation shows, with the load distri- 
bution and with the order of the harmonic. Deming (ref. 5) has shown by 
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calculation in special cases for the static propeller that the ring 
approximation is a reasonably good one. An effective radius of the order 
of 0. 8 r is considered reasonable for adoption in initial calc ula tions. 

The magnitude of the root -mean-square pressure is of interest, 

since most sound-recording instruments are calibrated in terms of it. The 
contribution of the harmonic of order m to the root -mean -square pres- 
sure is 


or 


rms 



+ 



■ rms 


= iL ( A 2 + B 2 )V2 

8* 2 


(26) 


where 


• 2*1 


Tx 


A = / ■{— cos(mB0 + ka) + 

J 0 


T ~(m + I) - Q 5® 

P 2 ' s / Rp2 


sin(mB0 + ka) 


d0 

S 


- 2 * 


B = 


Tx 

— sin(mB0 + ka) + 
S 2 




cos(mB9 + ka) 


d9 

s 


Wien rotational symmetry exists as for the condition of symmetrical 
loading, it is convenient in making the numerical calculations to let 
the field point be in the xy -plane, an arrangement which can always be 
attained by suitable choice of the line 0=0. 

Pressure relations- for the far field . - A further simplification in 
the results can be reached if the distance s from the propeller disk 
(0>yi^l) to the field point (x,y,0) is large. Then from 


s 


2 



\2 2 

- yi) + z i 


there is obtained 


S o “ s . 


yy± 


(27) 
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where s Q = V?+ y 2 (the distance from the field point to the center 
of the propeller disk). It also follows from 


S 2 = x 2 + p 2 (y - y d ) 2 + P 2 z d 2 


that 



where 

S D = \/x 2 + p2y 2 


( 28 ) 


If terms of order 1/s 2 are neglected in comparison with terms of 
order l/S, the sound pressure due to thrust is obtained from equa- 
tion (2k) by use of equations ( 5 ) and ( 28 ): 


icot 


0 — . . OC (S 0+ MX) n 2jt -imBe+ik-l-RpCOse 

P„~T^^e 3" i| (m + ^)J o e S 0 6 de ( 29 ) 


. k 

1 _1 «2 


From the known integral relation 

■ n2it 


0 


e i7\cos0-in9 de = 23 ti n J n (X) 


(30) 


it follows that 


tot -i^(S 0 tMt) 


±030 _ n 

p t ' T is; 1 e p 

Similarly, from equation ( 25 ), 


f( M + (51) 


P Q = 




R e 2 2* S o 


(52) 
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Observe that the argument of the Bessel function may be replaced by 


kyR e 

So 


- "“rot | - 
b o 


(33) 


where M rQ ^ ~ — the Mach number corresponding to the rotational speed 

at the effective radius. (The far-field approximation and the effective - 
radius approximation made use of in arriving at eqs. ( 31 ) and (32) need 
not, of course, be made simultaneously, since the required integration 
of the Bessel functions with respect to the radius may be carried out, 
at least numerically, without difficulty. ) 

Introduce the fundamental frequency denoted by = BO, so that 


a) = mBft = kc = 


The pressure magnitude for any harmonic m is then given by 


P = 


P T + P Q 


mcj^i 

2rtcS, 


t/m + - 

\ SJ B 2 


Q 


Be 


( kyR e 


J mB 


(34) 


This result may be compared with that for the case of zero forward speed 
(M = 0) given in reference 1: 



(35) 


It can be noted that the forward-speed Mach number affects the 
torque term containing Q only in the replacing of s Q = + y 2 (the 

distance from the field point to the propeller hub) by the smaller dis- 
tance S Q = This substitution occurs both in the argument 

of the Bessel function and in the outside factor l/S 0 . The thrust term 

containing T is more strongly affected, since, in addition to this 
change, the field effect ahead of and behind the propeller disk is influ- 
enced by the term M + ^ corresponding to a backward shift by ^- = -M 

increase because of the factor i/p 2 . 


and also by an 
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APPLICATION TO A SPECIFIC PROPELLER 


In order to give some indication of the effect of forward-speed 
Mach number on the sound pressure of a propeller, calculations based on 
equation (26) for the near field and on equation (34) for, the far field 
have been made. For this purpose, a two-blade research propeller having 
a 10-foot diameter and operating under the various conditions summarized 
in table I was chosen. As may be noted from this table, the propeller 
is assumed to operate at constant power j that is, the power coeffi- 
cient Cp and the torque Q are held constant. As the forward-speed 
Mach number is changed, the blade angle (3p> and consequently the 

thrust T are changed according to the propeller charts of reference 6 
so as to make the assumed conditions consistent with actual test oper- 
ating conditions. Only the fundamental, the first harmonic m = 1, is 
considered in these examples and the value chosen for the effective 
radius R e is 0.8R or 4 feet. It should be pointed out that the sound 
pressure computed from the data of table I would be obtained in pounds 
per square foot, since these data are given in English units. These 
values have been converted to dynes per square centimeter in the results 
given in the figures by multiplication by the approximate conversion 
factor 4 80 . 


Calculations for the near field . - Calculations of the root-mean- 
square pressure = -^^|p|^ based on equation (26) are made for 


various values of x in the range from — 0*5E to 0*5b along the line 
y = 0.6D (6 feet), that is, along a line l/5 of the propeller radius from 
the tip and extending a distance of 1 radius behind to 1 radius ahead of 
the plane of rotation. These results are shown plotted as a function of 
x/D for several Mach numbers in figure 4. It can be noted in this fig- 
ure that, for each Mach number, two peaks generally appear, one ahead of 
and one behind the plane of rotation. For the set of conditions under 
consideration (see table i), the highest peak amplitudes of pressure 
occur just behind the propeller plane at values of x/D in the range 
from -0.15 to -0.075. As the Mach number is increased from 0 to 0.4, 
the peak amplitudes decrease in magnitude, but as the Mach number is 
increased from 0.4 to 0.9 this trend is reversed and the peak amplitudes 
increase in magnitude, the peak for M = 0.9 being about 1.4 times that 
for M = 0. Also, as M increases from 0 to 0. 9 > the point at which 
the highest peak pressure occurs moves somewhat nearer the propeller 
plane. Thus, the generally severe sound -pressure conditions at take-off 
(H«0) tend to be alleviated in flight at the lower Mach numbers but may 
be reached again and even exceeded at the higher flight Mach numbers. 

The calculations for low Mach numbers presented in figure 4, it should 
be noted, are in substantial agreement with both calculated and measured 
results of reference 2. The high peak pressures obtained for the highest 
Mach numbers indicate that, with propeller -driven airplanes, the sound 



NACA TN 3018 


19 


pressures generated near the tips of the propellers are of significance 
and of possible concern with regard to both structural considerations 
and passenger comfort. 

Calculations for the far field . - Calculations of the root -mean - 
square pressures based on equation (3*J- ) for various values of x along 
the line y = 2D (20 feet) are shown plotted in figure 5* The trend with 
regard to Mach number in this figure is about the same as was noted in 
figure 4 for y = 0. 6 d, but the relative effect of the forward Mach num- 
ber appears greater at a distance. However, because of the greater dis- 
tances in figure 5 than in figure 4, the peak amplitudes are considerably 
less. 

Directional characteristics of the sound field . - For some purposes, 
especially with regard to calculations for the far field, it is desirable 
to consider the character of the sound pressure in terms of polar coordi- 
nates. For this purpose, substitutions can be made in equation (3*0 
(and eq. (26)) as follows: 

x = s Q cos $ y = s 0 sin $ s Q = ]Jx 2 + y 2 (36) 

The quantities x/S Q and y/S Q in equation (3*+ ) (and eq. (26)) 
may then be written as 



With these substitutions, equation (3*0 becomes 



(38) 

Calculations based on equation (38) for the Mach numbers of 0 and 
0.8 for a constant value of s Q = 2D are shown plotted as dashed curves 
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in figure 6. The solid curves in this figure represent results of calcu- 
lations along the line y = 2D obtained by replacing s n by — — — 

sin $ 

with y = 2D in equation (30) for the various Mach numbers considered 
in figures 4 and 5» A comparison of results of calculations along the 
line y = 2D with those along the circle s Q = 2D for M = 0 and 

M = 0.8 indicates that, for high Mach numbers (of the order of M = 0. 8), 
the peak pressures calculated along a line y = Constant are about the 
same as those calculated along the circle with radius equal to the con- 
stant value of y. Observe the second pressure peak which has developed 
in the forward location at M = 0.9* 

Separate components due to torque and to thrust of the sound field . - 
In order to give some indication as to the nature and proportion of 
sound pressure associated with each of the quantities, thrust T and 
torque Q, the root-mean-square pressures associated with each of these 
quantities for M = 0 and M = 0.8 are plotted for y = 0.6 d (6 feet) in 
figure 7. These plots show, as do equations (26) and (34), that the root- 
mean-square pressures associated with the propeller torque are symmetri- 
cally distributed with respect to the plane of the propeller for all Mach 
numbers, whereas the root-mean-square pressures associated with the 
thrust are symmetrical with respect to this plane only for M = 0. For 
the particular propeller and operating conditions under consideration, 
the amplitudes of pressure associated with thrust for low Mach numbers 
are higher than those associated with torque, but for hign Mach numbers, 
the opposite is true. In the interpretation of figure 7, it should be 
recalled that the results depend on the assumed operating conditions and 
that the torque Q is the same (2,680 ft-lb) in both parts of the figure; 
whereas the thrust is 1,850 pounds at M = 0 and 310 pounds at M = 0.8. 
Plots of this type can be used to obtain the sound pressure for various 
thrust and torque coefficients for a given propeller, since these coef- 
ficients appear as factors in equations (26) and (34) and hence can be 
normalized. 

It is of particular interest to note in the plot for M = 0. 8 that 
the peak pressures associated with both thrust and torque are consider- 
ably greater than the peak pressure associated with the sum of these 
quantities. This result indicates that the phase relationship between 
the two components is, in this case, such that each has a canceling 
effect on the other. Also in this case, the pressure associated with 
thrust may have its greatest value ahead of the plane of the propeller; 
whereas the total pressure has its greatest value behind the plane of 
the propeller. 

The discussion of the figures based on the specific examples illu- 
strates the fact that the results do depend markedly on the assumed 
operating conditions. Moreover, only the fundamental harmonic (m = l) 
for the specific two-blade propeller has been illustrated. Trend studies 
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on effects of higher harmonics, number of blades, and different oper- 
ating conditions would be of considerable interest. Some preliminary 
calculations on the sound pressures associated with the fundamental of 
a six-blade propeller have shown a greater relative effect of the 
forward-speed Mach number. 


CONCLUDING REMARKS 


Expressions have been given for the sound-pressure field due to the 
distribution of thrust and torque for any given harmonic of a rotating 
propeller in uniform subsonic flight. The general expressions (eqs. ( 17 ) 
to ( 23 )) involve integrations over the propeller disk; approximate 
expressions for the near field (eqs. (24) to ( 26 )) and for the far field 
(eqs. ( 31 ), ( 32 ), and (34)) involve integrations over a ring with an 
effective radius. The numerical examples have illustrated some free- 
space sound -pres sure results for the fundamental of a specific two -blade 
propeller under various operating conditions at various forward speeds. 

It is pertinent to remark again on some of the limitations of the 
analysis. The analysis presented utilizes the torque and the thrust 
distributions, which may be given empirically or theoretically, in such 
a way as to require that they arise purely from pressures acting normal 
to the surface of the blades. Empirical values of the thrust and torque 
include a contribution, generally small, due to skin friction and to 
separated flow. Some caution is then needed in the use of the results. 
For example, the empirical torque term is -somewhat larger than the torque 
due purely to normal pressures and, hence, the associated sound result 
due to torque may be overestimated slightly; similarly, the empirical 
thrust term may slightly underestimate the sound due to thrust. At high 
tip speeds there may be significant contributions to the noise due to 
wave drag associated with the thickness. ' These contributions are mainly 
taken into account by the effects of the wave drag on the torque. (Other 
sound effects, of the thickness should be separately calculated and 
included but these effects have not been explicitly presented herein. ) 

Another way of looking at the approximation, and perhaps a generally 
more convenient one for the study of trends, is to consider that the 
assumed torque and thrust distributions are actual theoretical ones, 
obtained if necessary by adjustment of the blades of the propeller in a 
potential flow; and hence the sound-pressure field is that corresponding 
to the chosen thrust and torque, but in potential flow. Although the 
sound-pressure field and the aerodynamic velocity field have been con- 
sidered as separate, it is of interest and of significance that the same 
concepts leading to the calculation of the sound field (made use of in 
the form of the acceleration potential) can lead to the linearized aero- 
dynamics of the propeller in compressible flow, including the represen- 
tation of the vorticity left behind in the wake as a result of spanwise 
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variation of loading. The theoretical induced drag and the theoretical 
wave drag (which is a form of acoustic loss) are inherently included in 
the representation employed. 

It may "be worthy of repetition that the pressure formulas as given 
tend to overestimate the contribution of the higher harmonics and that 
appropriate factors based on chordwise loading can be devised and may 
be required. Moreover, because of the theoretical (as well as the 
empirical) change in aerodynamic loading along the blade radially with 
increasing tip Mach number and forward-speed Mach number, the appropriate 
factors will change. In addition, the effective radius will be altered 
with the loading. Thus, although 0.8 blade radius may be suitable as an 
effective radius for normal loading, a smaller value for the effective 
radius may be more suitable for conditions where unloading of the tip 
occurs. 

A calculation study of trends under different assumed conditions, 
the effect of unsymmetrical loading, or of dual propellers, effects of 
obstacles or boundaries on the free -space results, and experimental 
confirmations for the in-flight propeller are interesting matters for 
further investigation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
' Langley Field, Va. , August 3^ 1955* 
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APPENDIX 

TEE SOUND-PRESSURE FIELD ASSOCIATED WITH 
A UNIFORMLY MOVING FORCE 
Derivation of Equations (4) and (6) 


Lamb (ref. 3, p. 502) derives from hydrodynamical principles the 
differential equation satisfied by the pressure condensation s* in the 
acoustic field associated with arbitrarily varying fixed forces acting 
on the medium 


#_ + al + 

dx 2 dy 2 dz 2 


2 at 2 / 


= jJM + Sr + dz\ 

c 2 \3x dy dz / 


(Al) 


where Iamb's X, Y, Z are actually the extraneous forces per unit volume 
divided by the density p. (in effect. Lamb is dealing here with the 

acceleration potential. ) In terms of the perturbation pressure p = pc 2 s* 
the equation is 


d^p = div F 



dx 3y dz 


(A2) 


where F is the arbitrary force per unit volume having components 
F„, F, r , F„. Lamb shows that if the periodic force Ae ia)t is imagined 

x y z / \ • 

concentrated on an infinitely small space at ( x i>yi> z i) and to be in the 

direction of x, the pressure at the field point (x,y,z) (the distance s 
from the location of the force ) is given by 


P * 



(A3) 


where 


s = (j(x - x i ) 2 + (y - yi ) 2 + ( z - z i ) 2 
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so that the concentrated force is equivalent to a double source or 
acceleration doublet whose axis is in the direction of the force. For 
a general harmonically periodic force F having components F x , Fy, F z 

varying in time as e icut 


P * 



V) 


e lks 

s 





(A4) 


where k = ^. 

c 

For an arbitrary time dependent concentrated force F(t) located 
at the pressure at field point (x,y,z) can be expressed as 


P = 


4it s 


1 

a F *(‘ - !) 

, a - 1) 

, a -t) 


dx s " 

8 y s 

c)z s 


(A5) 


where the differentiations affect the result only through the variable s. 

The extension to a moving concentrated force may, within the frame- 
work of small-perturbation theory, be made in several ways - for example, 
as mentioned in reference 7 . One formal procedure utilized by Kiissner, 
for example, makes use of the invariance properties of the wave equation 
and utilizes combined Galilean and Lorentz transformations. The following 
direct procedure given briefly in reference 7 is believed to be of 
intrinsic interest, as it represents a simplification in Prandtl's pro- 
cedure (ref. 8 ) for the case of moving constant source distributions, a 
procedure which consists of scheduling a succession of fixed sources in 
a path to act consecutively one after the other so as to represent in 
effect the desired source moving along its path. let the arbitrary con- 
centrated force act only as an impulse during an infinitesimal interval 
at time t = r. The impulse may be written as 


F(t) 6 (t - t) 


(a6) 
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where the impulse function &(t) = 0 for t ^ 0 and is characterized 
for t = 0 as having unit area with respect to t. The useful property 
of the impulse function of "sifting or selecting" a value of a function 
is exhibited by the following relation (see, for example, ref. 9, p. 6l): 



t ) 5 (T)dT 


f(r)&(cr - t )dT 


- t (o) 


(A7) 


Let a succession of such impulses act, one following the other, in a path, 
points along which are given in space -fixed coordinates by 

I = I(t) tj = t](t) g = £(t) 

The effect at time t of all such impulses which act before the time t 
is then given from Lamb's result (eq. (A5)) as 



where, as has been defined, there is a nonzero contribution to the 
integral only for values of t defined by the characteristic relation 

t - T ■ | = |||[x - + [y - ^I] 2 + 5 - £(0] 2 (A9) 

which expresses the distance between the source point and the field point 
in terms of the time of travel of the outgoing waves. 

The integral in equation (A8) corresponds to a summation of temporary 
fixed sources. To represent the case of uniform rectilinear motion with 
velocity V in the positive x-direction, let the sources be located on 
the |-axis and flow consecutively one after the other at the positions 

I = Vt n = 0 C=0 (A10) 

so that for t = -« the source is located at | = -» and for t = 0 
the source is located at the origin. It follows that the distance 
between source point and field point is 


s = \|(x - Vt)^ + y^ + z^ 


(All) , 
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Replace in equation (A8) the variable r by 0 where 


t - T 



(A12) 


By the sifting property of the 5-function (eq. (A7 ) ) the integral will 
have a value only for 9=0 and equation (A8) becomes 


P = 


idivBiisi 

bn |_ 8 


0=0 ' 


(A13) 


where the quantity within the bracket is now to be evaluated. Equa- 
tions (All) and (A12) give a quadratic equation for t 'which, on choice 
of the solution that leads to t < t and with 0=0, gives 



where 

S = {{x - Vt) 2 + p 2 (y 2 + z 2 ) 

0 = 'll - Vp 



c 


(Al4) 


Equations (All) and (A12) also give 


1 dT _ c 

s d0 cs _ v(x - Vt) 


where replacing s by its characteristic value c(t - t) from equa- 
tion (A9) and replacing r by its value given in equation (Al4) yields, 
for 9=0, 


1 dT _ 1 

s d9 ’ S 


(A15) 


where S is defined in equation (Al4)l Equation (A13) becomes 


P = 





S 


(Al6) 
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This result is referred to space-fixed coordinates; it is finally 
desired also to express this result in terms of a field point ( x o^0^ z o) 

of a coordinate system moving uniformly along with the source located at 
x l> Vi* Z 1 so that one may put 


x - vt = x 0 - x-l y = y 0 - yp z = z o - z i 


After this substitution is made, the zero subscripts may be dropped to 
yield as the end result corresponding to equation (A5) for the source 
and the coordinate system in uniform rectilinear motion: 


P = 



(ait) 


where 

s = ij(x - x-l) 2 + 3 2 [(y - yi) 2 + (z - zp )*\ 

M(x - x^) + S 


p = \/l - M 2 


Equation (A17 ) corresponds to the result given in equation (4) of 

the analysis. For the periodic harmonic force F varying as it 

may he expressed in the form 


i — 0 -i ka 

p = i ( F • V) — (A18) 

4jt S 

which is equivalent to equation (6) of the analysis. Equations (AI 7 ) 
and (Al8) are the sought-for generalizations of equations (A5) and (A4) 
of the appendix^ for the case where the disturbance and the field point 
are in uniform rectilinear motion. Comparison of these equations shows 
at once that s is replaced by S in the amplitude and by a in the 
phase. A geometric interpretation of these quantities is shown in fig- 
ure 1 and discussed in the analysis following equation ( 5 ). 
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Remarks On the Case of the Moving Disturbance 
and the Fixed Observer 

A few remarks are' in order on the significance of equation (Al 6 ). 

In this equation, the field point (x,y,z) is given in space-fixed coordi- 
nates while the disturbance force (or propeller) is in the in-flight 
condition moving with velocity V in the positive x-direction. Its 
location is given by x^_ = Vt, y^ = 0, = 0 so that at t = 0 it is 

located at the origin. It may be observed that the distance between 
disturbance and field point S is numerically the same whether given 
by the relation in equation (Al 6 ) or (AI 7 ). Hence, the pressure magni- 
tude at any observer location for a sound -radiating element of the 
in-flight propeller is the same as that for the observer moving along 
with the propeller, provided the proper instantaneous distance between 
observer and propeller is used. However, there will be a difference in 
the frequency perceived by the observer. This frequency will be that 
for the case of the uniformly moving observer modified by the Doppler 
effect} thus, the frequency of each harmonic is modified by c/(c + V r ), 
where V r is the component of the propeller -forward speed (with proper 

sign: minus, for approaching, plus for receding from, the observer) 
in the acoustical direction from the observer to, say, for the far 
approximation, the hub of the propeller. The acoustical direction is 
not quite that from the observer to the propeller location; it actually 
points from the observer to the location of the propeller when the sound 
which reaches the observer was emitted (direction QP rather than Q0 
in fig. l). The Doppler frequency factor c /( c + V r ) is given geomet- 

QU Cj 

rically in figure 1 by the ratio of — = — . 

QP cr 
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Figure 1.- Geometric representation of s, S, and a 
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(a) Sharp pulse. 

(b) Triangular pulse. — — = 0.03. 

2rtr 

(c) Rectangular pulse. = 0.03.- 

2itr 

(d) Rectangular pulse. — = 0.06. 

2otr 

Figure 3*- Effect of impulse shape on the relative amplitudes 

of the harmonics. 
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Figure 7-- Distribution of the root-mean-square total pressures and the root-mean-square pres- 
sures associated with thrust and torque for a two-blade, 10-foot-diameter propeller at forward- 
speed Mach numbers of 0 and 0.8. Cp = 0.10; m = 1; k = 0.29686; y = 6 feet. (Operating con- 
ditions of propeller given in table I.) 



